CONSTRUCTING ALL MAGIC SQUARES OF ORDER THREE 



GUOCE XIN 

Abstract. We find by applying MacMaiion's partition analysis that all magic squares of 
order three, up to rotations and reflections, are of two types, each generated by three basis 
elements. A combinatorial proof of this fact is given. 
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1. INTRODUCTION 

A magic square of order n is an n by n matrix with distinct nonnegative integer entries 
such that every row sum, column sum, and (two) diagonal sums equals to the same number m, 
the magic number. Adding 1 to every entry will give us a traditional magic square of positive 
integers. A magic square is pure if the entries are the consecutive numbers from to — 1, 
and hence it has magic number 3 ("3 ^) • ^ 

Magic squares have been objects of study for centuries. As Pickover wrote in his book^ p. 
60]: 

. . . the holy grail of magic squares creation would be to discover a method that would 
generate every possible arrangement for a square of a given size. Such a solution is 
probably not discoverable. 

This "holy grail" could be achieved by first finding the complete generating function (which is 
a rational function) for magic squares of a given size, and then writing the generating function 
as a sum of simple rational functions, the series expansion of which has only nonnegative 
coefficients. 

We achieve this for magic squares of order 3, as given in Theorem 12.11 
Weak maqic squares, magic squares without the restriction of distinct elements, have been 
studied in by using the rich theory of counting solutions of a system of linear 

Diophantine equations, or equivalently, counting lattice points of a convex polytope. For 
further references, see p, Ch. 4.6]. These methods also apply to counting magic squares, but 
give no obvious reason why a simple solution as in Theorem 12. II exists. 
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2 GUOCE XIN 

We give our main result in Section 2, and give a combinatorial proof in Section 3. In Section 
4, we discuss the discovery of our main result and possible future work. 



2. Main Results 

A magic square of order 3 is a 3 by 3 matrix of distinct nonnegative integers such that every 
row sum, column sum, and diagonal sum equals the magic number m. 

Our main result is the following Theorem \'2.1\ which generates all magic square of order 3. 
Let 
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(2.1) 



Then they are related as follows: 



B = C + D, Ti=B + C, T2 = B + D. (2.3) 

If we let C be obtained from C by reflecting in the vertical axis, then we have one more relation: 
D = C + C. It is straightforward to check that A, C, and D are linearly independent. 

In fact, A,C,D are the three basis elements that generate all magic squares of order 3, 
and Ti,T2 are the unique magic squares with magic numbers 12 and 15, respectively, up to 
rotations and reflections. 

Theorem 2.1. Every magic square of order three, up to rotation and reflection, can be written 
uniquely as either Ti + iA + jB + kC or T2 + iA + jB + kD, where i,j, k are nonnegative 
integers and A, B,C, D,Ti,T2 are as in ()2.H1 . ()2.2j) . 

Remark 2.2. Note that traditional magic squares can be generated by either iA + jB + kC 
or iA + jB + kD for positive integers i,j, k. This description reveals a kind of symmetry. 

Theorem 12.11 savs that magic squares, as a set of lattice points, is a disjoint union of 16 = 8-2 
polyhedrons that are isomorphic to N'^, where the factor 8 is the order of the dihedral group of 
rotations and reflections. We will give a combinatorial proof of this result in the next section. 
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Corollary 2.3. The number of magic squares of order 3 with magie number 3s and its asso- 
ciated generating function is given by 
8t\l + 2t) 



:i-t)(l-t2)(i-t3) 



J2 (25^ - + 1 - (-1)' + 3)) 

s>0 ^ ^ 

= 8 (i^ + 3t^ + + 7t^ + lOt^ + + 17i^° + 22t^^ + 26t^^ + ■•■). 

3. A Combinatorial Proof 

In what follows, magic squares are always of order 3 unless specified otherwise. 
Let M be a magic square with magic number m. We write 
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(3.1) 
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where 

Cl: Every row sum, column sum, and diagonal sum is equal to m. 
C2: The entries of M are distinct nonnegative integers. 

Rotating or reflecting M will give us different magic squares. Without loss of generality, we 
can assume that C3 is smaller than ai, 03 and ci, and that ci < 03. Also by subtracting A times 
the minimal entry of M from M, we can assume that is an entry of M. Then M satisfies 
the following two extra conditions: 

C3: One of the entries of M is 0. 
C4: C3 < ai,a3,Ci, and Ci < a^. 

In fact, C4 can be replaced with 

C4': C3 < Cl < 03 < oi, 

which follows from the sum of the two diagonals. 

If M satisfies the above four conditions, then we say that M is a reduced magic square. It 
is well-known that the magic number m is m = 362- Let m = 3s or equivalently s = 62- 

Lemma 3.1. If M is a reduced magic square, then 02 = 0, and C2 — 2s. 

Proof. Since b2 = s = m/3, where m is the magic number, the last statement C2 = 2s follows 
from a2 = 0, which is what we are going to show now. 

In a reduced magic square M, ai and are the largest two entries among the four corners 
fli) ^3) Cl, C3. It follows from the first and third row sums and column sums that 02 < &i, C2, 63. 
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We see that &2 = s > 4, since all entries are distinct nonnegative integers. It remains to 
show that C3 cannot be 0. Assuming that C3 equals 0, then ai = 2s by the diagonal sum 
fli + &2 + C3 = 3s. By investigating the first row sum and the first column sum, we get 
as < s — 1 and ci < s — 1, contradicting the condition for the diagonal (03, 62, ci). □ 

Lemma 3.2. A reduced magic square M can be uniquely written as Ti + aC + (3D, where 
a > — 1 and /3 > are integers. 

Proof. To see the existence, we use Lemma l3.ll Assuming that C3 = r and 62 = s, we obtain 
all the entries of M by the condition CI for row sums, column sums, and diagonal sums: 



M 



2s — r s + r 

2r s 2s — 2r 
s — r 2s r 



Comparing the above matrix with 



Tl + aC + (3D 



7 + 2a + 3/? 5 + a + 3/5 

2 + 2(3 A + a + 2(3 6 + 2a + 2/3 
3 + a + /? 8 + 2a + 4/3 1 + /3 



we solve uniquely for a and (3: 



a = s — 2r — 2, and (3 = r — 1. 



Consequently, 



s = a + 2/3 + 4 and r = (3 + 1. 



We see that C3 > 1 and Ci > C3 implies that a > — 1 and /3 > 0, completing the proof of the 
existence. 

The uniqueness follows from the above proof, and also from the fact that C and D are 
linearly independent. □ 



We are now ready to give the proof of our main theorem. 
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Proof of Theorem, \2.1\ It is straightforward to check that 
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Ti + iA+ jB + kC 



T2 + iA + jB + kD 



7 + i + 5j + 2k i 5 + i + Aj + k 

2 + 1 + 2] A + i + d>i + k Q + i + Aj + 2k 

?, + i + 2j + k % + i + Qj + 2k l + i+j 

8 + z + 5j + 3A; i 7 + ii + 4j + 3fc 
A + i + 2j + 2k 5 + i + 3j + 2A; Q + i + Aj + 2k 

'i + i + 2i + k 10 + 2 + 6j + 4A; 2 + i+i + k 



(3.2) 



(3.3) 



give different magic squares for all nonnegative integers j, k. 

Given a magic square M, we need to show that M equals either (j3.2p or ()3.3p . 

Let i be the minimum of the entries of M. Then up to rotations and reflections, we can 
assume M' = M — iA is a, reduced magic square. By Lemma f3. 21 M' can be uniquely written 
as Ti + aC + l3D, with a > -1 and /? > 0. 

li a > p > 0, M' can be rewritten (recall that B = C + D) as Ti+ pB + (a - p)C. Hence 
we let j = /3 > and = a - /3 > 0. 

If a < (3, M' can be rewritten (recall that Ti + D = T2 + C) as 

Ti + aB + {p- a)D = T2 + C + aB + {[3 - a - l)D = T2 + {a + l)B + {(3 - a - 2)D. 

Thus we let j = a + 1 > and k = (3 — a — 2 > —1. 

The only remaining case is A; = — 1, which is equivalent to /3 = a + 1. But in this case 



M' = Ti + {[3- l)C + I3D 



5 + 5/3 4 + 4/3 
2 + 2/3 3 + 3/3 4 + 4/3 
2 + 2/3 6 + 6/3 1 + /3 



which is not a magic square because it has equal entries. 



□ 



4. Further Discussion 

The combinatorial proof in the previous section seems unlikely to be applicable to magic 
squares of higher order. We describe how we discovered Theorem 12.11 by using MacMahon's 
partition analysis, which has been restudied by Andrews and his coauthors in a series of papers 
(see e.g., 

MacMahon's idea is to use new variables to replace linear constraints. For example, if we 
want to count nonnegative integral solutions of the linear equation ai + a2 — = 0, we can 
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simply write the generating function as 

^ ^ ^ ^ A 12 3^ 1-Axi 1-AX2 I-X3/A ' 

ai,a2,a3>0 ai,a2,a3>0 ^ ^' ' 

where CTa means to take the constant term in A. Then the counting problem is converted to 

evaluating the constant term of a special rational function, which can be done be computer 
n 

package as in For a rigorous description about how the above works in general situation, 
i.e., in a field of iterated Laurent series, the reader is referred to 0]. 

Using a computer we can easily obtain the generating function of weak magic squares of 
order 3: 

(1 — tXA:X,'jX'iX(,X2Xj,X^X^X\) (1 + tXi,X-jX<iX^X2XzX^X%X\f' 



G 



[1 — tXiX^XgX^'^Xs'^Xs'^) (1 — ^72:5X3X42X2^X92) 

1 



X 



[1 — tXyXsXsXi^Xg^Xg^) (1 — tXiX5X9X72x22x62) , 

where t records m/3 since the m is always divisible by 3, and the exponents in xi, ...Xg 
represents ai, 02, 03, 61, ... . 

To obtain the generating function for magic squares, we shall take only terms in G that have 
different exponents in the x's. To eliminate those terms with same exponents in xi and X2, we 
subtract by the diagonal diag^.^ ^^G with respect to Xi and X2, where 



rGN sGN rGN 

and we use the formula for a rational power series F{x, y): 

diag.,F(x,y) = CT ^--i^F(A,,A2). 

Similarly, we can eliminate those terms with same exponents in Xj and xj for all i and j. 

The generating function of all magic squares of order 3 is still complicated. We can add the 
extra constraints that C3 < ci < 03 < ai to eliminate rotations and reflections. It suffices to 
find a way to add the constraint that the exponent of Xg is smaller than that of X7. The other 
constraints can be added iteratively. We omit the details here. 

Finally we obtain the generating function of desired magic squares: 

t^X7^X5'^X3^Xi ''X8^X6^X9X42 (1 + tXiX5X9X42x8^X32 — 2t2x52a;9X42x8'^Xi^X3^X7X6^) 



tX7X5X3Xi2x82x62) (1 — tX4X7X9X6X2X3X5X8Xi) 

1 

< 

(1 — t'^X^'^XgXi'^Xs'^Xi^Xs^XYXQ'^) (1 — t'^X72x53xi^X8^X3'^X6^X9X42 



(4.1) 
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We observe that part of the numerator can be rewritten as 
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1 I J. 2 2 2 r,j.2 2 2 4 3 3 2 

1 + 1X1X5X9X4 Xg X3 — 2t X5 X9X4 Xs Xi X3 X7X6 

/I j2 2 2433 2\ij. 222/1 j 222\ 

= (1 — r X5 X9X4 xg xi X3 X7X6 ) + 1x1X5X9X4 Xs X3 (1 — 1x7X5X3X1 xg xe ) 



where both terms on the right-hand side will cancel with the denominator of ()4.1|1 . Theorem 
12.11 then follows. 

The order 4 case would be really hard. The difficulty lies in the fact that there are 880 pure 
magic square of order 4 (up to rotations and reflections), which suggests that there will be 



at least 880 simple rational functions. Our current package as provided in 







is not powerful 



enough to find an explicit generating function for magic squares of order 4 analogous to ()4.1|) . 
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